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t-H 1 Abstract. We present new results of 3-D AMR MHD simulations focusing on two distinct as- 
pects of PPN evolution. We first report new simulations of collimated outflows driven entirely by 
magnetic fields. These Poynting flux dominated "magnetic towers" hold promise for explaining 
key properties of PPN flows. Our simulations address magnetic tower evolution and stability. 
We also present results of a campaign of simulations to explore the development of accretion 
disks formed via wind capture. Our result focus on the limits of disk formation and the range 

f^) ] of disk properties. 
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1. Introduction 

2 1 Jets are observed in proto-Planetary Nebuale (PPN), Young Stellar Objects, radio 

galaxies and other astrophysical objects. Models suggest that jets are launched and 
collimated by accretion, rotation and magnetic mechanisms (Pudritz et al. 2007). One 
\ paradigm to explain the formation of jets in PPN is that of a binary system in which 
i—i ■ material from an AGB star is accreted onto an intermediate-mass companion (Soker & 
' Rappaport 2000; Nordhaus, Blackman & Frank 2007). Jets or collimated outflow have 
(N. recently been observed in the close binaries NGC 6326 and NGC 6778 (Miszalski et al. 

2011). Previous numerical studies of wind capture in binaries have found enhancements 
^-j- ■ over Bondi-Hoyle accretion rates onto the secondary (Mastrodemos & Morris 1998; Pod- 
siadlowski & Mohamed 2007; De Val-Borro et al. 2009). These results, however, must 
be tested using high resolution simulations because the implications for the maximum 
\ outflow power are dramatic; the answer could rule in, or out, the secondary as the engine 
i— I ■ powering jets in PPN. 

In the case of magnetized jets, plots of specific angular momentum vs. jet poloidal 
speed are distinct for different MHD engines (Ferreira et al. 2006). This fact can help 
constrain the physics of the jet engine. The importance of the magnetic fields relative to 
the flows' kinetic energy, divides jets into (i) Poynting flux dominated (PFD; Shibata & 
Uchida 1986), in which magnetic fields dominate the jet structure, (ii) magnetocentrifugal 
(Blandford & Payne 1982), in which magnetic fields only dominate out to the Alfven 
radius. The observable differences between PFD and magnetocentrifugal jets are unclear, 
as are the effects of cooling and rotation on PFD jets. Recently, cooling magnetized jets 
have been formed in laboratory experiments (Lebedev et al. 2005). Such studies, along 
with high resolution 3-D MHD numerical simulations, have proven to play a key role in 
the understanding of the physics of jet launch and stability. 
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2. Magnetic towers 

2.1. Model 

We use the Adaptive Mesh Refinement (AMR) code AstroBEAR2.0 (Cunningham et al. 
2009; Carroll-Nellenback et al. 2011) to solve the equations of radiative-MHD in 3D. 
The grid represents 160 x 160x400 AU divided into 64x64x80 cells plus 2 AMR levels. 
Initially, the molecular gas is static and has an ideal gas equation of state (7 =5/3), 
a number density of 100 cm s -1 and a temperature of 10000 K. The magnetic field is 
helical, centrally localized and given by the vector potential (in cylindrical coordinates) 
A(r, z) = [(r/4)(cos(2r) + l)(cos(2z) + l)]</> + [5 (cos(2r) + l)(cos(2z) + l)]k, for r, z < 
30 AU, and A(r, z) = elsewhere. The magnetic pressure is higher than the thermal one 
inside the magnetized region. 

Source terms are used to continually inject magnetic or kinetic energy at cells r, z < 
30 AU. We carry out 4 simulations: an adiabatic, a cooling (Dalgarno & McCray 1972) 
and a rotating (Keplerian) PDF jet, as well as a hydrodynamical jet which has the same 
time average propagation speed and energy flux than the adiabatic PFD jet. 

log(n) [cm -3 ] 
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Figure 1. Logarithmic false color number density (cm -3 ) maps of the gas at the midplane of 
the cubic computational domain. Adiabatic (1st column), rotating (2nd column) and cooling 
(3rd column) PDF jets. The hydrodynamic jet is in the 4th column. From top to bottom the 
time is 42, 84 and 118 yr. 



2.2. Results 

Magnetic pressure pushes field lines and plasma up, forming magnetic cavities with low 
density. The adiabatic case is the most stable. PFD jets decelerate relative to the hydro- 
dynamic one; magnetic energy pressure produces axial but also radial expansion. PFD 
jet cores are thin and unstable, whereas the hydrodynamic jet beam is thicker, smoother 
and stable. The PFD jets are sub-Alfvenic. Their cores are confined by magnetic hoop 
stress, while their surrounding cavities are collimated by external thermal pressure. PDF 
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jets carry high axial currents which return along their outer contact discontinuity. Pinch, 
m = 0, kink, m = 1, and m = 2 current-driven perturbations are evident in the PDF jets. 
The perturbations are amplified by cooling, firstly, and base rotation, secondly. This hap- 
pens because both shocks and thermal pressure support are weakened by cooling, and 
the total pressure balance at the jets' base is affected by rotation. See Huarte-Espinosa 
et al. 2011 (in prep) for details. 

3. Binary-formed disks 

3.1. Model 

We use Astrobear2.0 (Cunningham et al. 2009; Carroll-Nellenback et al. 2011) to solve 
the equations of hydrodynamics in 3D. The grid represents 100 AU 3 divided into 32 3 cells. 
We also use 5 AMR levels which allow us to have high resolution, of order 0.4 AU, around 
the model stars. These, are implemented with gravity particles, separated by 25 AU 
and follow circular orbits. The primary simulates an AGB star with a mass of 1.5 M©, 
and a spherical constant wind of lOkms" 1 and M ~ 10~ 8 M© yr _1 . The secondary star 
simulates a main sequence star or a white dwarf with 1 M© . 




-3.0-2.0-1.0 0.0 1.0 2.0 3.0 
^ 10 AU 

Time=2,9 (orbit) 

Figure 2. Logarithmic number density (cm -3 ) map of the gas in the orbital plane. The primary 
and secondary stars are located at (0.6,-0.9) and (-0.9,1.2), respectively. The blue and white 
components are the AGB wind, the spiral structure of which is caused by the orbital motion. 
The red gas component is the accretion disk that is formed around the secondary via wind 
capture. 

Initially, we implement an isothermal (7 = 1.001) gas in the grid which is unaffected 
by the gravitational field of either of the stars, but they pull each other around the 
center of gravity which is located at the origin (Figure 2). The system orbits twice and 
the grid is filled with material from the primary's AGB wind. The orbital motion causes 
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a helical-like distribution on the wind. We then "turn-on" the effect of the secondary's 
gravity on the gas and allow the system to orbit ten times under these conditions. 

3.2. Results 

As the stars orbit each other, the dense gas from the AGB wind flowing near the secondary 
is captured by the star's gravitational field. This process is essentially a form of Bondi- 
Hoyle flow and the flow pattern is localized near the secondary. A spiral shock pattern 
does however propagate through the AGB wind. 

Within one orbital timescale we see that an accretion disk forms around the secondary. 
We note that previous 3D simulations of accretion disk formation via wind capture in 
binaries have explored star separations smaller than our model with 25 AU (e.g. Mas- 
trodemos & Morris 1998). Thus our results show that disks can form even out to large 
separations. 

Preliminary analysis finds the density and temperature of the disk progressively in- 
creasing in time. We also find that disk material close to the secondary follows Keplerian 
orbits confirming it is bound to the star. See Huarte-Espinosa et al. 2012 (in prep) for 
details. 

4. Conclusions and Summary 

PFD jet beams are lighter, slower and less stable than kinetic-energy dominated ones. 
We predict characteristic emission distributions for each of these. Current-driven pertur- 
bations in PFD jets are amplified by cooling, firstly, and base rotation, secondly: shocks 
and thermal pressure support are weakened by cooling. Total pressure balance at the 
jets' base is affected by rotation. Our simulations are in good agreement with the models 
and experiments of Shibata & Uchida 1986 and Lebedev et al. 2005, respectively. 

Regarding the problem of the accretion disk formed by wind capture, we see the for- 
mation of accretion disks at large binary separations of order 25 AU. The density and 
temperature of the disk progressively increases in time. 
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